Production and decay of Sulphur excited species in a ECRIS plasma 
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The most important processes for the creation of S 12+ to S 14+ ions excited states from the ground 
configurations of S 9+ to S 14+ ions in an electron cyclotron resonance ion source, leading to the 
emission of K X-ray lines, are studied. Theoretical values for inner-shell excitation and ionization 
cross sections, including double KL and triple KLL ionization, transition probabilities and energies 
for the deexcitation processes, are calculated in the framework of the multi-configuration Dirac- 
Fock method. With reasonable assumptions about the electron energy distribution, a theoretical 
Ka X-ray spectrum is obtained, which is compared to recent experimental data. 
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I. INTRODUCTION 

Electron-cyclotron resonance ion sources (ECRIS) are 
characterized by their capacity to produce large popu- 
lations of highly charged ions and by high electron tem- 
peratures. X-ray emission, including brcmsstrahlung and 
characteristic lines, caused by these high-energy electrons 
due to electron-ion collisions, has been used for plasma 
diagnostics. 

In 2000, Douysset et al. [l[ proposed to estimate the 
ionic density of each charge state in an ECRIS plasma 
through the measurement of the intensity of the emitted 
Ka lines. As pointed out by these authors, this calcula- 
tion depends, on the following conditions: 

• precise identification of the X-ray lines emitted by 
the different charge states; 

• a deep understanding of the different paths leading 
to the excited states and their cross sections; 

• knowledge of the shape of the electronic distribu- 
tion function above the threshold for production of 
excited states. 

In 2001, some of us published an analysis of K X- 
ray spectra emitted by Ar ions in an ECRIS plasma [l|, 
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showing that a complete analysis of these spectra calls 
for a careful examination of all excitation and ionization 
processes that lead to the excited states of the different 
ionic species whose decay will yield the detected lines. We 
showed that singlc-ionization processes were not enough 
to explain all the observed features and only when dou- 
ble KL ionization is taken into account is it possible to 
account for all of them. 

Several ECR ion sources with superconducting coils 
and superconducting or permanent magnet hexapoles are 
now becoming available 0, 0j H, @| ■ One of them, the 
Electron-Cyclotron resonance ion trap (ECRIT) at the 
Paul Scherrer Institute [1] , has been specifically designed 
for X-ray spectroscopy of the ions in the plasma. This 
ECRIT has a large mirror ratio, a large plasma vol- 
ume, and is an ideal instrument for detailed observa- 
tion of the plasma for highly-charged medium- Z ions. 
Even though the frequency of the injected microwave 
was low (6.4 GHz), X-rays of core-excited ions ranging 
from heliumlike to quasi neutral sulfur, chlorine and ar- 
gon have been measured. Thanks to this source and a 
high-resolution spherically-bent crystal X-ray spectrom- 
eter high-resolution spectra (« 0.3 eV) are now available 
0, HI, which represent an order of magnitude improve- 
ment over the previous experiment |l[ This improved 
resolution calls for a theoretical study more refined than 
the one we used for the previous work, and allow for the 
identification of new mechanism, as more lines can be 
resolved. 

In the present work, we examine the more important 
atomic processes that may take place in the PSI ECRIT, 
and contribute to the creation of low-lying excited states 
with a K hole which will de-excite by the emission of K 
X-ray lines, in order to find out the ionic density of each 
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charge state in a plasma. We used a multi-configuration 
Dirac-Fock (MCDF) program, including intra- and outer- 
shell correlation, relativistic and quantum electrodynam- 
ics (QED) contributions, to obtain values of energies and 
transition probabilities for the subsequent de-excitation 
processes. Theoretical values for excitation and ioniza- 
tion cross sections, including double and triple ionization 
in K- and L-shells, were computed. 

Combining these results with estimated values of ion 
densities, we were able to generate theoretical Ka X-ray 
spectra. A comparison with an experimental spectrum 
has led to more realistic values of ion densities. 

This paper is organized as follows. In Sec. HU we 
present briefly the MCDF method, with emphasis on the 
specific features of the code which are important for en- 
ergy correlation calculations. In Sec. IIII) we analyse 
sulphur spectra obtained in an ECRIS detailing the ex- 
periment, the peak energies calculations, the processes 
leading to the creation of ionic excited states, and the 
calculation of the corresponding cross sections, and fi- 
nally the calculation of the line intensities. In Section llVl 
we present the results and discuss them. 



II. THE MCDF METHOD 

The general relativistic MCDF code developed by J. 
P. Dcsclaux and P. Indclicato 0, [LH was used to 
calculate bound-state wave functions and energies. De- 
tails of the method, including the Hamiltonian and the 
processes used to build the wave functions can be found 
elsewhere 

The total wave function is calculated with the help of 
the variational principle. The total energy of the atomic 
system is the eigenvalue of the equation 

M n ° pait ^n,J,u(- •• .ri,. ..) = -Bn.j,M*n../.A/(- ■ ■ ,r ? , . . .), 

(1) 

where II is the parity, J is the total angular momentum 
eigenvalue, and M is the eigenvalue of its projection on 
the z axis J z . In this equation, the hamiltonian is given 
by 
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where TCd is the one electron Dirac operator and Vij is an 
operator representing the electron-electron interaction of 
order one in a. The expression of Vij in Coulomb gauge, 
and in atomic units, is 
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where ry = |r,- — r ^ | is the inter-electronic distance, u>ij 
is the energy of the exchanged photon between the two 
electrons, at are the Dirac matrices and c is the speed of 
light. We use the Coulomb gauge as it has been demon- 
strated that it provides energies free from spurious con- 
tributions at the ladder approximation level and must be 
used in many-body atomic structure calculations [Til . [lH ] . 

The term Q represents the Coulomb interaction, the 
term ((4]) is the Gaunt (magnetic) interaction, and the 
last two terms (0 stand for the retardation operator. In 
this expression the V operators act only on r,j and not 
on the following wave functions. 

The MCDF method is defined by the particular choice 
of a trial function to solve the Dirac equation as a linear 
combination of configuration state functions (CSF): 
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The CSF are also eigenfunctions of the parity 77, the total 
angular momentum J 2 and its projection J z . The label v 
stands for all other numbers (principal quantum number, 
etc.) necessary to define unambiguously the CSF. The 
c v are called the mixing coefficients and are obtained by 
diagonalization of the Hamiltonian matrix coming from 
the minimization of the energy in Eq. ([1]) with respect to 
the c v . 

The CSF are antisymmetric products of one-electron 
wave functions expressed as linear combination of Slater 
determinants of Dirac 4-spinors 
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where the ip-s are the one-electron wave functions and 
the coefficients di are determined by requiring that the 
CSF is an eigenstate of J 2 and J z . The di coefficients are 
obtained by requiring that the CSF are eigenstates of J 2 
and J z . 

The Multi-Configuration approach is characterized by 
the fact that a small number of configurations can ac- 
count for a large amount of correlation. 

The so-called Optimized Level (OL) method was used 
to determine the wave function and energy for each state 
involved. In this way, spin-orbitals in the initial and final 
states for the radiative transitions are not orthogonal, 
since they have been optimized separately. This non- 
orthogonality effect is fully taken into account [l6|, H3] , 
using the formalism proposed by Lowdin [lj| . The length 
gauge has been used for all radiative transition probabil- 
ities. 

Radiationless transition probabilities were calculated 
using Desclaux's code [HI]. The bound wave functions 
were generated using this code for configurations that 
contain one initial inner-shell vacancy while the contin- 
uum wave functions were obtained by solving the DF 
equations with the same atomic potential of the initial 



3 



state. With this treatment, the continuum wave func- 
tions are made orthogonal to the initial bound state wave 
functions, thus assuring orthogonality. The continuum 
wave function is normalized to represent one electron per 
unit time. 

QED contributions for the transition energies and tran- 
sition probabilities have been found to be negligible. For 
example, the QED contributions for the ls2s — > Is 2 tran- 
sition energy and probability are, respectively, 0.03% and 
0.01%. 



III. ANALYSIS OF SULPHUR SPECTRA IN AN 
ECRIS PLASMA 

A. Experiment 

Sulphur X-ray spectra were obtained in the 
Electron 6.4 GHz Cyclotron Resonance Ion Trap 
(ECRIT) at the Paul Scherrer Institut (PSI) [if, 
by the Pionic Hydrogen Collaboration (see 
http : //www . f z- j uelich . de/ikp/ exotic- atoms/ index 
The experimental set-up was composed mainly of two 
parts: the ECR ion trap (the X-ray source) [4 
Bragg spectrometer set-up in Johann geometry 
The X rays reflected by the spectrometer crystal were 
recorded by a two-dimensional position-sensitive detec- 
tor [23[ placed in the proximity of the Rowland circle of 
the spectrometer. 

In order to obtain X-ray spectra from highly-charged 
sulphur, SO2 gas was injected. To improve the ionization 
efficiency, a gas mixture with O2 was used, adjusted to 
a mixing ratio of about 1:9 (the main gas is oxygen). 
Details about calibration and spectra construction can 
be found in Refs. [H HI- 

The obtained spectra cover an energy region corre- 
sponding to He-like, Li-like, and Be-like Ka line energies 
of sulphur ions. 



ing the high-precision QED calculation from Ref. [26[ . In 
this case the theoretical results differ from the experimen- 
tal ones by 0.46 eV and 0.41 eV for the two Li-like lines, 
respectively, as shown in Tab. [I] and in Fig.[TJ These dis- 
crepancies are due to the higher accuracy of the present 
experimental data. Moreover the improved resolution 
(0.3 eV) in comparison with the Ar case (3 eV) makes the 
comparison between theoretical and experimental spec- 
tra more sensitive to shifts. Therefore a new calculation 
was performed, including correlation up to the 4/ sub- 
shell. For that purpose, the basis space was enlarged to 
include now all orbitals in the n = 3, 4 electronic shells. 
A much better agreement with experiment is obtained 
(see Tab. HJ, as the energy differences between theory 
and experiment were reduced to 0.06 eV and 0.04 eV, re- 
spectively, for the two Li-like lines, thus showing the im- 
portance of electronic correlation in the interpretation of 
high resolution experiments. We arc dealing here with 
auto-ionizing initial states in these transitions, except in 
the case of heliumlike line. An auto-ionizing state corre- 
spond to a resonance degenerate with a continuum, cor- 
pkjjipponding to a final state with a free electron, which 
lead to both a broadening (due to the Auger transition 
probability) and a shift. Up to now this shift has been 
evaluated only for neutral atoms with an inner-shell hole 
[27l . [28j . Here, a rough estimate of this effect shows it 
must be in the order of a few meV to a few tenth of eV, 
depending on the number of allowed Auger channels. 

Using our calculation, including only correlation up to 
n = 4, we obtain for the Ml line the energy 2430.28 eV, 
0.07 eV below the experimental value. A more accurate 
energy shift could be deduced, if the experiment was cali- 
brated against our MCDF value, as some cancellations in 
the correlation energy should occur: both the initial and 
final states of the main Li-like transitions, for example, 
correspond to the Ml line initial and final states with an 
extra electron. 



B. Transition energies 

The spectra obtained at PSI cover the 2.400-2.460 keV 
energy range. The more important features in these spec- 
tra are the Be-like ls2s 2 2p 1 P 1 -> ls 2 2s 21 S El line 
at 2.418 keV, the He-like ls2s 3 Si -> ls 21 S Ml line 
at 2.430 keV, and the Li-like \s2s2p 2 Pj -> ls 2 2s 2 S 1/2 , 
J = 1/2,3/2 lines, at 2.437 keV and 2.438 keV, respec- 
tively. 

In a preliminary calculation, we have used the MCDF 
code, including only intra-shell correlation, to calculate 
the energies of these lines in the same way as in our pre- 
vious calculations for Ar ions spectra Q . The basis space 
includes all possible electron configurations built from Is, 
2s, and 2p orbitals corresponding to single and double ex- 
citations from the main configuration. The experimental 
spectrum energy scale has been fixed by setting the He- 
like ls2s 3 Si -> Is 2 ^0 Ml transition to 2430.351 eV, us- 



C. Excited species production 

Radiative spectra obtained in ECRIS systems usually 
include radiative lines resulting from de-excitation of ions 
in different charge states and different levels. A detailed 
analysis of the main processes leading to those initial 
states is needed. We refer the reader to Martins et al. @ 
for details. 

We assume that all S q+ ions, where q is the degree of 
ionization (q = Z — m, m being the number of electrons 
in the ion), are initially in the ground configuration. The 
two main processes leading to a hole in the K shell are 
K ionization of the S^ 9_1 ^ + ion ground configuration and 
excitation of the S 9+ ion ground configuration. In order 
to generate a preliminary theoretical spectrum, we took 
into account processes that create single K-hole excited 
states with two, three, and four electrons from the ground 
configurations of sulphur ions. 

The intensity of the line corresponding to a transition 
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from level i to a level j in a S 9+ ion with a K hole and 
in level i, is given by 

/-,ir r vr". (8) 

In this equation fkv and A q - are the radiative i j tran- 
sition energy and probability, respectively, and N^ ,q is 
the density of S 9+ ions with a K hole and in level i. This 
density is obtained from the balance equation 

N q (NeV^—^ + N*- 1 (N^af^- 1 ^) = Nf' q A q . 

(9) 

Here, N e is the electron density, Nq and Nq are the 
densities of S 9+ and S^ 9-1 -^ ions, respectively, in the 
ground configuration, A q is the transition probability for 
de-excitation of these ions by all possible processes (ra- 
diative and radiationless), a f~ CKC ' q [ s the excitation cross 
section for the processes leading from the S 9+ ion in the 
ground configuration to the excited level i of the same 
ion with a K hole, and of Ion ^ 9 i s the ionization 
cross section leading from the S^ 9-1 ^ ion in the ground 
configuration to the excited level i of the S q+ ion with a 
K hole. 

The quantities {N e v e Oi} give the rate of the number 
of events related to a process (excitation or ionization) , 
averaged over the electron distribution energy, and are 
defined by 

/>00 

(N e v e ai)=N e f e (E)v e (E)<j(E)dE, (10) 
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where f e (E) is the electron distribution energy, and 
v e (E) and o(E) are, respectively, the electron velocity, 
and the cross section for the process, for electron energy 
E. 

Electron-impact excitation cross sections were com- 
puted using the first Born approximation following the 
work of Kim et al. 0, [3(| • In these calculations we ob- 
tained the cross sections for the processes leading from 
each level j of the S 9+ ion ground configuration, to the 
excited level i of the S 9+ ion with a K-shell hole, Oji. We 
used a Multi-Configuration Dirac-Fock (MCDF) wave 
function for the atom and a Dirac wave function for the 
free electron. Only the Coulomb interaction between the 
free electron and the atomic electrons was considered. 
The individual cross sections thus obtained were then 
weighted by the statistical weight gj of each j level of 
the ground configuration in order to obtain the excita- 
tion cross section o i 

K-shell ionization cross-sections were computed using 
the relativistic binary-encounter-dipole (RBED) model 
[3ll ]. This method has provided accurate results [Hj], 
even at energies close to the ionization threshold. We use 
a simplified version of this method, referred as the rela- 
tivistic binary-encounter-Bethe (RBEB) model, in which, 
besides the incident electron kinetic energy T, only three 
orbital constants are needed: the kinetic energy of the 
ejected electron U, the orbital binding energy B, and the 
electron occupation number N for the pertinent shell. 
The RBEB cross section expressions reads 
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where t = T/B, u = U /B, v t is the speed of an electron 
with kinetic energy T, Vb is the speed of an electron with 
kinetic energy B, v u is the speed of an electron with ki- 
netic energy U, a is the fine-structure constant, ag (0.529 
A) is the Bohr radius and 

Pt = v t /c, # = l-_L-_ t' = T/mc\ (12) 
(3 b = v b /c, P 2 b =l - 7jjyj2 > b' = B/mc 2 , (13) 
Pu=v u /c, 01 = 1- * , u' = U/mc 2 . (14) 

(i + u y 

The product of the ionization cross section, leading 
from the S'- 9-1 -'"*" ion in the ground configuration to the 



S 9+ ion with a K hole, by the statistical weight gi of 
the level i, yields the K-shell ionization cross section 
of lon,( - 9 l ' q ' leading from the S^ 9-1 ^ ion in the ground 
configuration to the excited level i of the S 9+ ion with a K 
hole. We assumed a statistical population for each level 
(with a particular J value) in the final configuration. 

As the electron energy distribution is not known, we 
used Eq. ijSJ) and @, with cross sections computed by 
us. We assumed an electron energy distribution which is 
a linear combination of a Maxwellian distribution (90%), 
describing the cold electrons at the thermodynamical 
temperature (fcT=lkeV), and a non-Maxwellian distri- 
bution (10%), describing the hot electrons (fcT=20keV) 
heated by electron cyclotron resonance j33[. Details of 
the calculation will be given elsewhere. The ion densities 
provided by Douysset et al. [l| for the corresponding 
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Ar ions were used as starting values. We note that the 
initial level of the He-like Ml transition originates from 
the Is 2 He-like ions ground configuration by excitation 
and from the ls 2 2s Li-like ions ground configuration by 
single ionization. The initial levels of the two El Li-like 
transitions also originate from the same configuration by 
excitation. 

Although the theoretical spectrum thus obtained spans 
the 2378-2460 eV energy range, it is shown in Fig. [1] in 
the 2420-2455 eV energy range for comparison with ex- 
perimental data (dash-dotted line.) The calculated spec- 
trum is normalized to the 2438 eV peak intensity. As 
Fig. [T] shows, the theoretical spectrum resulting from ex- 
cited states obtained from the ion ground states by ex- 
citation and single K-ionization processes only is enough 
to account for the major features in the experimental 
spectrum. In Tab. |TT] we list the ion densities used in 
the theoretical spectrum, obtained through an iterative 
adjustment to the experimental spectrum. 

From the comparison with experiment we note that, 
although the Li-likc line intensities are well reproduced, 
the Ml He-like line theoretical intensity falls below the 
experimental intensity. At the same time, some minor 
features remain unexplained. So, we looked for processes 
that could lead to other excited states suceptible of ex- 



plaining the discrepancies between our calculated values 
and the experimental data. 

We found out that contributions from the following 
configurations had also to be included: S 13+ ls2p 2 , and 
S 12+ \s2s2p 2 , and ls2p 3 . We took into account all exci- 
tation, single-ionization, double KL-ionization and triple 
KLL-ionization processes from the ground state of S 9+ to 
S 14+ ions, leading to these configurations. A diagram of 
all excitation, ionization and decay processes considered 
in this work is shown in Fig. [2] 

As an example, the S 12+ ls2s 2 2p 1 P 1 excited configu- 
ration can be obtained through Is — > 2p excitation from 
the S 12+ Is 2 2s 2 configuration, through K ionization from 
the S 11+ \s 2 2s 2 2p configuration, through KL double- 
ionization from the S 10+ ls 2 2s 2 2p 2 configuration, and 
through KLL triple- ionization from the S 9+ \s 2 2s 2 2p i 
configuration. 

In what concerns the calculation of the double KL and 
triple KLL ionization cross-sections, we used the semi- 
empirical formula developed by Shevelko and Tawara 
[34| . with the fitting parameters proposed by Belcngcr 
et al. [35| . and the ionization energies calculated in this 
work (sec Tab.|lH|. 

Taking in account all the referred processes, Eq. ([8]) 
now reads 
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In this equation, erf 1 ""™' and <7f LL ~ lon,( " 9 3 ' q ^ are the 

double and triple ionization cross sections respectively 



IV. RESULTS AND DISCUSSION 

The theoretical spectrum, including all lines resulting 
from the decay of the excited states shown in Fig. [21 ob- 
tained using the methods discussed above and assum- 
ing, for each line, a linear combination of a gaussian 
and a lorentzian distributions, designed to approximate 
a Voigt profile with 0.3 eV width, is presented in Fig. [3] 
in a semilogarithmic scalc(a), and in a linear scale(b). In 
Fig. [3^a) is shown, as a dash-dotted line, the theoretical 
spectrum resulting from excited states obtained from the 
ion ground states by excitation and single K-ionization 
processes only. This is enough to account for the ma- 
jor features in the experimental spectrum, clearly seen in 
Fig. [3](b) . Therefore we looked for other processes that 
could lead to new excited states also from the ground 
states. Fig. |3£a) shows that inclusion of doublc-KL and 



triple-KLL ionization allows for a much better theoretical 
interpretation of these features. 

If only single ionization and excitation processes are 
considered, the ls 2 2s and Is 2 ground configurations are 
the only ones that contribute to the experimental sulphur 
spectrum in the 2420-2455 eV energy range, the first one 
through single ionization, and the second one through 
monopolc excitation. 

In this case, the 2446.68 eV and 2448.28 eV lines, re- 
sulting from the decay of the ls2p 3 Pi 2 levels are not 
significant, as these levels are poorly fed by single exci- 
tation from the Is 2 ground configuration. Also based in 
the single excitation and ionization processes, the calcu- 
lated intensity of the Ml line resulting from de-excitation 
of the ls2s 3 Si level is lower than the corresponding ex- 
perimental line intensity, after adjusting the calculated 
spectrum to the 2438 eV line intensity. 

However, the ls2s configuration may also be obtained 
through double- and triplc-ionization processes, from the 
ls 2 2s 2 and ls 2 2s 2 2p configuration states, respectively 
Including the contribution of these processes in the cal- 
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culatcd spectra allows for a much better fit to the exper- 
imental spectra. 

The intensity ratio between the Mf ls2s 3 Si — > Is 2 1 S 



J 



and the ls2s2p 2 P 3 / 2 — > ls 2 2s 2 Si/ 2 lines (denoted i and 
2 in Fig. [T]) is given by 
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x g(ls2s 3 5'i). (16) 



Here N q+ arc the densities of S 9+ ion in the ground con- 
figuration and g is the statistical weight of the ls2s 3 Si 
level, and TYi is the transition yield of line i. 

The Is2s2p 2 ion configuration is obtained when the 
double KL and triple KLL ionization processes are taken 
in account. The ls2s2p 23 S 1 -> ls 2 2s2p 3 Pj, J = 0,1,2 
triplet visible at the left in the figure in now predicted in 
agreement with experiment. 

One should note that there are several weak features in 
the spectrum in Fig. [3j which are not reproduced. These 
features, are slightly above statistical noise. Some of 
them could be due to n = 3 or 4 satellite lines of the 
hcliumlikc and lithiumlikc ions. Yet our calculations did 
not produce any line that could explain the observed fea- 
tures. A search in X-ray transition energy database [28| 
could not produce any line that could explain those fea- 
tures. 

The energy and transition yield values calculated in 
this work for all lines contributing to the theoretical spec- 
trum are listed in Tab. HVl The excitation cross sections 
are shown in Table [V] for 10 keV impact energy, as an 
example. 



V. CONCLUSIONS 

Ion densities in the plasma depend on the electron 
energy distribution. For hot electrons this distribution 
is non-Maxwcllian, making difficult the calculation of 
the ion densities. So, for the generation of the theo- 
retical spectrum, we started from the ion densities pro- 
vided by Douysset et al [l|. Individual line intensities 
were obtained using Eq. (JSJ) with the cross section val- 
ues computed by us assuming a linear combination of 
Maxwellian distribution, describing the cold electrons, 



and a non-Maxwellian distribution, describing the hot 
electrons, and the estimated ion densities. After compar- 
ison with the experimental spectrum, ion densities were 
then adjusted in order to match the experimental peak 
intensities. In this way we were able to obtain a reason- 
able agreement between theory and experiment. 

In this work we have shown evidence of the contribu- 
tion of several new mechanism that contribute to the cre- 
ation of excited state in an ECRIS plasma. In particular 
double and triple processes must be considered to repro- 
duce properly line intensities. We have also shown that 
MCDF transition energies with medium size configura- 
tion space can reproduce reasonably well most features 
of the spectra, including energies and branching ratios. 
In some part of the spectrum there are hint of weak lines 
that could not be reproduced so far, despite extensive 
search, including n = 3 satellites of heliumlike and lithi- 
umlike lines. 
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TABLE I: Electronic correlation effect on the Is 2s 3 Si — > Is 2 
So Ml line energy and on the energy shifts, relative to the 
Ml line, of the Is 2s 2p 2 P 1/2 -> Is 2 2s El(l) and Is 2s 2p 
2 P3/2 — * Is 2 2s El(2) lines, compared to experiment. The 
energy values are in eV. 





Ml 


El (1) 


El (2) 


Intra-shell correlation (up to 2p) 


2430.09 


6.30 


7.34 


Outer-shell correlation (up to 4f) 


2430.29 


6.82 


7.79 


Experiment 




6.76 


7.75 



TABLE II: Ion density values obtained through an iterative 
adjustment to the experimental spectrum. 

Ion densities (m -3 ) 

7TZ 



gl0+ 
gll+ 
g!2+ 
g!3+ 
g!4+ 



8.0x J 
7.0X 1 
7.0X 1 
9.0X 1 
1.4X 1 
7.0X 1 



TABLE III: Calculated energies for single K- (Is 



KL- (ls^ 1 2s- 1 
ization (in eV). 



double 

2\ 



Is" 1 2P" 1 ) and triple KLL- (Is" 1 2s~ 2 ) ion 



glO+ 

S n+ 

gl2+ 
gl3+ 
g!4+ 



Configuration 



ls" 



2s 2 2p 3 
2s 2 2p 2 
2s 2 2p 
2s 2 
2s 



Is" 1 2s" 1 Is" 1 2P" 1 Is" 1 2s" 



2799.63 
2880.79 
2967.42 
3059.32 
3135.92 
3222.42 



3332.02 
3476.37 
3626.01 
3787.89 
3929.29 



3328.31 
3472.45 
3619.21 



4344.02 
4608.55 
4869.29 
5141.14 



TABLE IV: Energy and Transition yield (TY) values cal- 
culated in this work for all lines contributing to the the- 
oretical spectrum. 



q 


conf; 


LSJi 


conff 


LSJ f 


E (eV) 


TY 


13+ 


Is 


2 

2s" 


^Si jo 


i 2 

Is 


2p 


2 P , 

r 3 /a 


2378.45 


1.75E-02 


13+ 


Is 




Si /o 


Is 


2p 


ll /2 


2380.36 


1.02E-02 


12+ 


Is 


2p 3 


5 S 2 


Is 


2p 2 


3 Pa 


2387.45 


1.06E-03 


12+ 


Is 


2p 3 


5 s 2 


Is 


2p 2 


3 Pi 


2388.45 


7.10E-04 


12+ 


Is 


2p 3 


3 D 3 


Is 


rn 2 

2p 2 


'D 2 


2397.23 


1.97E-03 


12+ 


Is 


2p 3 


3 D 3 


, 2 

Is" 


„ 9 

2p 2 


3 P 2 


2404.12 


1.27E-01 


12+ 


Is 


2p 3 


'Pi 


^ 2 

Is 


rn 2 

2p 2 


'So 


2404.20 


3.39E-01 


12+ 


Is 


2p 3 


3 D 2 


-i 2 

Is 


n 2 

2p 2 


3 P 2 


2404.39 


7.82E-02 


12+ 


Is 


2p 3 


3 Dr 


Is 


rn 2 

2p 2 


3 P 2 


2404.39 


8.66E-03 


12+ 


Is 


2s 2p 2 


"Da 


-i 2 

Is 


2s 2p 


'Pi 


2405.31 


3.93E-01 


12+ 


Is 


r* 3 

2p 3 


3 D 2 


Is 


r» 2 

2p 2 


3 Pi 


2405.40 


3.80E-01 


12+ 


Is 


2p 3 


3 Dr 


-i 2 

Is 


2p 2 


3 Pi 


2405.40 


1.29E-01 


12+ 


Is 


2s 2 2p 


3 Pi 


, 2 

Is 


2s" 


'So 


2405.75 


3.37E-03 


12+ 


Is 


2s 2p 2 


3 Pi 


i 2 

Is 


2s 2p 


'Pi 


2406.00 


3.82E-02 


12+ 


Is 


r* 3 

2p 3 


3 D! 


, 2 

Is 


r» 2 

2p 2 


3 Po 


2406.04 


2.31E-01 


12+ 


Is 


2s 2p 2 


3 P 2 


i 2 

Is 


2s 2p 


'Pi 


2407.05 


1.93E-01 


12+ 


Is 


r* 3 

2p 3 


'D 2 


, 2 

Is 


r» 2 

2p 2 


'Da 


2408.18 


7.78E-01 


12+ 


Is 


2p 3 


3 Si 


i 2 

Is- 


2p 2 


3 P 2 


2409.99 


3.73E-01 


12+ 


Is 


2p 3 


3 P 2 


, 2 

Is 


2p 2 


'Da 


2410.23 


2.99E-01 


12+ 


Is 


2p 3 


3 Si 


, 2 

is- 


n 2 

2p 


3 Pi 


2411.00 


2.82E-01 


12+ 


Is 


2s 2p 2 


3 D 3 


, 2 

Is- 


2s 2p 


3 Pa 


2411.10 


3.28E-01 


12+ 


Is 


2s 2p 2 


3 D 2 


i 2 

Is- 


2s 2p 


3 Pi 


2411.41 


2.94E-01 


12+ 


Is 


2p 3 


3 Si 


, 2 

Is 


n 2 

2p 2 


3 Po 


2411.64 


9.44E-02 


12+ 


Is 


2s 2p 2 


3 Pi 


i 2 

Is 


2s 2p 


3 Pa 


2412.07 


2.75E-01 


12+ 


Is 


2s 2p 2 


3 D 2 


, 2 

Is 


2s 2p 


3 Pa 


2412.22 


9.74E-03 


12+ 


Is 


2s 2p 2 


3 Po 


i 2 

Is 


2s 2p 


3 Pi 


2412.29 


5.29E-01 


12+ 


Is 


2s 2p 2 


3 Dr 


Is- 


2s 2p 


3 P 2 


2412.39 


1.39E-01 


12+ 


Is 


2s 2p 2 


3 P 2 


Is- 


2s 2p 


3 P 2 


2413.04 


8.62E-01 


12+ 


Is 


2s 2p 2 


3 Pi 


^ 9 

Is- 


2s 2p 


3 Pi 


2413.27 


3.83E-01 


12+ 


Is 


2s 2p 2 


3 Di 


i 2 

Is 


2s 2p 


3 Pi 


2413.36 


2.70E-01 


12+ 


Is 


2s 2p 2 


3 P 2 


, 2 

is- 


2s 2p 


3 Pi 


2414.23 


2.76E-02 


12+ 


Is 


2s 2p 2 


3 Pi 


_ 9 

Is- 


2s 2p 


3 Po 


2414.24 


3.69E-02 


12+ 


Is 


2p 3 


"Da 


, 2 

Is" 


r» 2 

2p 2 


3 P 2 


2415.07 


1.75E-01 


13+ 


Is 


2s 2p 


4 Pl/2 


^ 9 

Is 


2s 


2 Si/2 

/ 


2415.18 


5.28E-01 


13+ 


Is 


2s 2p 


4 P 3 /2 


i 2 

Is- 


2s 


Si ii 


2415.69 


2.14E-02 


12+ 


Is 


2s 2p 2 


'Pi 


- 9 

Is- 


2s 2p 


'Pi 


2416.02 


9.80E-01 


12+ 


Is 


2p 3 


'D 2 


- 9 

Is- 


2p 2 


3 P1 


2416.07 


1.66E-02 


12+ 


Is 


2p 3 


3 Pi 


- 9 

Is- 


„ 9 

2p 2 


3 P 2 


2416.49 


5.44E-01 


13+ 


Is 


2s 2p 


4 P 5 /2 


- 9 

Is- 


2s 


2 Si /2 


2416.93 


1.21E-01 


13+ 


Is 


2p 2 


4p 

r 5 j2 


- 9 

Is- 


2p 


2 P 3 /2 


2417.01 


9.90E-01 


12+ 


Is 


2p 3 


3 P 2 


., 9 

Is- 


2p 2 


3 P 2 


2417.13 


1.06E-01 


13+ 


Is 


2p 2 


4p 

/ 


- 9 

Is- 


2p 


2 Pl/2 


2417.16 


9.29E-03 


12+ 


Is 


2p 3 


3 Pi 


- 9 

Is- 


2p 2 


3 Pl 


2417.49 


1.26E-01 


12+ 


Is 


2p 3 


3 Po 


- 9 

Is- 


„ 9 

2p 2 


3 Pl 


2417.56 


1.00E+00 


12+ 


Is 


2p 3 


3 P 2 


- 9 

Is- 


2p 2 


3 Pl 


2418.13 


3.24E-02 


12+ 


Is 


2p 3 


3 P1 


- 9 

Is 


2p 2 


3 Po 


2418.14 


1.46E-01 


12+ 


Is 


2s 2p 2 


'So 


Is 2 


2s 2p 


'Pi 


2418.49 


4.64E-01 


12+ 


Is 


2s 2 2p 


'Pi 


Is 2 


2s 2 


'So 


2418.58 


3.41E-01 


12+ 


Is 


2p 3 


'Pi 


Is 2 


2p 2 


'D 2 


2420.98 


3.85E-01 


12+ 


Is 


2s 2p 2 


3 Sr 


Is 2 


2s 2p 


3 Pa 


2423.78 


3.72E-01 


12+ 


Is 


2s 2p 2 


3 Sr 


Is 2 


2s 2p 


3 Pi 


2424.97 


1.80E-01 


12+ 


Is 


2s 2p 2 


3 Sr 


Is 2 


2s 2p 


3 Po 


2425.94 


5.06E-02 


12+ 


Is 


2s 2p 2 


'Da 


Is 2 


2s 2p 


3 P 2 


2427.61 


1.50E-02 


12+ 


Is 


2s 2p 2 


3 Pi 


Is 2 


2s 2p 


3 P 2 


2428.29 


6.35E-01 


12+ 


Is 


2s 2p 2 


3 Po 


Is 2 


2s 2p 


3 Pi 


2428.68 


8.25E-02 


Continued on '. 


Next Page. 
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TABLE IV - Continued 



q 


confi 


LSJi 


conf/ 


L5J/ 


£ (eV) 


TY 


12+ 


Is 


2s 2p 2 


3 P 2 


Is 2 


2s 2p 


3 P 2 


2429.35 


1.19E-01 


12+ 


Is 


2s 2p 2 


3 Pi 


Is 2 


2s 2p 


3 Pi 


2429.48 


2.85E-02 


14+ 


Is 


2s 


3 Sx 


Is 2 




'So 


2430.29 


1.00E+00 


12+ 


Is 


2s 2p 2 


3 Pi 


Is 2 


2s 2p 


3 Po 


2430.45 


6.82E-02 


12+ 


Is 


2s 2p 2 


3 P 2 


Is 2 


2s 2p 


3 Pi 


2430.54 


2.88E-02 


13+ 


Is 


2p 2 


2 D 5/2 


Is 2 


2p 


2 P 3 /2 


2431.02 


3.53E-01 


13+ 


Is 


2p 2 


2 D3/2 


Is 2 


2p 


2 Pl/2 


2432.09 


9.94E-01 


13+ 


Is 


2p 2 


2 Pl/2 


Is 2 


2p 

— v 


P 3 /2 


2433.28 


2.80E-01 


13+ 


Is 


2p 2 


Pl/2 


Is 2 


2t> 

— F 


2 Pl/2 


2435.19 


6.57E-01 


13+ 


Is 


2p 2 


o ' 
P3/2 


Is 2 


2d 


P 3 /2 


2435.49 


9.22E-01 


13+ 


Is 


2s 2p 


2 Pl/2 


Is 2 


2s 


2 Si /o 

/ * 


2437.11 


9.39E-01 


13+ 


Is 


2p 2 


2 P 3 /2 


Is 2 


2p 


2 Pl/2 


2437.40 


7.83E-02 


13+ 


Is 


2s 2p 


2 P 3 /2 


Is 2 


2s 


2 Si/2 


2438.08 


1.00E+00 


14+ 


Is 


2p 


3 Pl 


Is 2 




'S 


2446.68 


1.00E+00 


13+ 


Is 


2s 2p 


2 P 1/2 


Is 2 


2s 


2 Si/2 


2447.24 


3.24E-01 


13+ 


Is 


2s 2p 


2 P 3 /2 


Is 2 


2s 


2 Si/2 


2447.64 


9.06E-01 


14+ 


Is 


2p 


2 3 P 2 


Is 2 




'So 


2448.28 


3.10E-01 


13+ 


Is 


2p 2 


Sl/2 


Is 2 


2p 


2 P 3 /2 


2449.62 


6.60E-01 


13+ 


Is 


2p 2 


Sl/2 


Is 2 


2p 


2 Pl/2 


2451.53 


1.89E-01 


14+ 


Is 


2p 


J Pl 


Is 2 




'So 


2460.15 


1.00E+00 



TABLE V: Sulphur ions excitation cross sections, for 10 keV 
electron impact energy. 



q 


initial conf. 


final conf. 


a(m 2 ) 




12+ 


Is 2 2s 2 'So 


Is 2s 2 2p d Pi 


2.05 x 10" 


-!27 




1 J2 r,„ 2o 


'Pi 

i _ n'A 2n 


3.05 x 10" 


-28 
■HI 



5l/2 



Is 2s 2p 



3l/2 

Pl/2 
2p(2) 

r i/a 

2p(l) 
r l/2 

4 P 3 /2 
2p(l) 
r 3/2 
2n(2) 

3/2 



6.68 x 10~ 29 
1.62 x 10 -25 
1.64 x 10~ 26 
3.68 x 10~ 28 
3.45 x 10" 25 
1.06 x 10~ 26 



14+ 


Is 2 


'So 


Is 2p 


3 Pi 


2.29 x 10-*' 










'Pi 


3.21 x 10" 25 










3 Si 


3.72 x 10~ 26 
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FIG. 1: Theoretical spectra for K excitation and single K 
ionization processes obtained with intra-shell ( ) and outer- 
shell correlation (- ■ -) plotted against experimental data (•). 
The peak 1 refers to the Ml Is 2s 3 Si — > Is 2 line, the peaks 
2 and 3 have contributions from the Is 2s 2p 2 Pi/2,3/2 - * Is 2 
2s and the Is 2p 3 Pi — > Is 2 lines, and the peak 4 refers to the 
Is 2s 2p 2 3 Si -> Is 2 2s 2p 3 Po,i,2 lines. 
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(9+) 1s 2 2s 2 2p 3 



■ ► 



Radiative transition 

Two electron-one photon radiative transition 
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K ionization 
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\ • 
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■i \ 
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. \ 
\ \ 
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v--^-vV--->Hs2s 2 2p 



1s2p 3 



> 1s 2 2p 2 
-► 1s 2 2s2p 
-► 1s 2 2s 2 



(13+) 1s 2 2s 




* 1s 2p' 
s 2s 2p 



1s 2 2p 
1s 2 2s 



(14+) 1s 2 




FIG. 2: Excitation, ionization and decay processes considered 

in this work: K excitation ( ); single K ionization (- ■ -); 

double KL ionization (- • • -); triple KLL ionization(- • • 
• -); radiative decay ( — ), and radiationless decay ( • ■ ■ ). 
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FIG. 3: Theoretical spectra obtained with single ionization 
(- • -) and multiple ionization ( — ) processes plotted against 
experimental data (•) in log scale (a) and linear sclae (b). For 
the legend of the peaks, see Fig. [T] caption. 



